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The complex potential of two cylinders moving in an infinite liquid is determined by the 
method of image doublets, and the solution is expressed as an infinite series in rectangular 
coordinates. Approximate solutions in finite form are given for various cases. A method for 
generalizing the solution for the case of more than two cylinders is indicated. Applications 
to the flow induced by a cylinder moving in the presence of plane boundaries are given and 
the stream lines are illustrated in certain cases. 



1. Introduction 

The problem of the motion of two cylinders in any 
manner with their axes always parallel was appar- 
ently first solved by Hicks. 1 The velocity potentials 
are found in general as definite integrals, which, 
when the cylinders move as a rigid body, are ex- 
pressed in finite form as elliptic functions of bipolar 
coordinates. Greenhill 2 has expressed the solution 
as an infinite series in bipolar coordinates. 

In a problem in wave motion the need arose for a 
solution of the potential flow for a cylinder moving 
parallel to a plane boundary. The problem is the 
same as that of two parallel cylinders of equal radii 
moving with equal velocities normal to the line of 
centers. Thus the flow is a special case of the classic 
results obtained by Hicks and Greenhill. However, 
it was necessary to have the potential function 
represented in rectangular coordinates to determine 
the additional term representing the wave motion. 
The bipolar coordinates used by Hicks and Greenhill 
are defined geometrically, and the transformation to 
rectangular coordinates does not admit a convenient 
representation. Therefore, the solution was ob- 
tained by the method of image doublets and is here 
generalized to the case of two cylinders moving in any 
manner with their axes always parallel. The present 
method may also be applied to the motion of more 
than two cylinders which was not done in the bipolar 
system with either Hicks' or GreenhhTs analysis. 

Only ideal fluids without vorticity are considered 
in this study. 

2. Solution for Two Cylinders 

Consider two parallel cylinders C and C of radii 
b and b' momentarily at a distance / apart moving 
normal to their axes in a liquid at rest at infinity. 
Choosing the coordinate system as in figure 1, the 
cylinder velocities may be expressed as Ue ia and 
U'c ia ' , respectively. 

The complex potential w = < t>o+i*Po due to C, in 
the absence of C, is that of a doublet whose axis is in 
the direction a 



w =Ub 2 — > 

z 



a) 



"This problem was encountered in connection with a wave Investigation spon- 
sored by the Office of Naval Research, and the present account is published in 
order to make more convenient formulas available. 

i W. M. Hicks, On the motion of two cylinders in a fluid, Quart. J. Pure and 
Appl. Math. 16, 113 (1879). 

2 A. G. Greenhill, Functional images in Cartesians, Quart. J. Pure and Appl. 
Math. 18, 231 (1882). 




Figure 1. Notation diagram for motion of two cylinders in an 
infinite Liquid. 



where z is the complex coordinate position. The 
normal velocities introduced on C by w are cancelled 
by the image doublet whose axis is in the direction 



where 



Wi= - Ub iil^(rf7) 

(by 



/.= 



/ 



(2) 
(3) 



The sum Wo-jrWi is the well known expression for the 
potential of a doublet in the presence of a cylinder 
(see Milne-Thomson 3 ). Similarly, the normal ve- 
locities introduced on C by w x are cancelled by the 
second image doublet 



where 



/.= 



6 2 



/-/. 



(5) 



Continuing, an infinite sequence of image doublets 
of decreasing strength is obtained. The general 
expression for w k is 



3 L. M. Milne-Thomson, Theoretical Hydrodynamics, Article 8.82, p. 221 
(MacMillan & Co., Ltd., London, 1938). 
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Wk=* 



e~ 



j. 



for k even 

a 

for & odd 



where p =Ub 2 , fc 



*-(/-.» 

-*,(*) 



(6) 



'\2 



(66') 2 



and 



* (/-/*- 2 ) 2 (/-/*-.) 2 
6 2 



ft- 2 



(7) 



/r 



f~~Jk-2 

{b'Y 



for fc even 
for fc odd 



(7a) 



it being understood / =0. 



Similarly, starting with the potential of C in the 
absence of C. 



Wo '=U'(b')*- 



z-f 



«-&<&($££ 



and 



pia' 



W k =■ 



**-(S-fi 



tt for k even 



-ia' 



w 



7 for A: odd 



(8) 
(9) 

(10) 



where 



^=u'(br,^=-u 



w@)' 



ft= 



(66') 2 



and 



(60 2 



/; 



*=-< 



6 2 



for A; even 



for & odd 



(11) 



(Ha) 



and again /o'^O- 

Thus the complex potential w=<j>-\-i\p for the 
motion of two cylinders is 



ff?=S (w*+w*) 



fc=0 



(12) 



where the w fc are obtained from eq (1) through (7) 
and the w k from eq (8) through (11). 



3. Approximate Solutions 

The series (12) converges rapidly for small values 
of bb ' lj 2 so that for this case one or two terms of the 
series gives a close approximate solution 



W 

w ^w +Wq +w 1 -\-w 1 f for -sr small. 

f 



(13) 



For larger values of W If 2 a more accurate solution 
may be obtained by approximating the remaining 
terms. For this purpose let 



?7=lim/ fc (k odd)=lim// (k even) 
7/^lim/fc (&even)=lim// (k odd) 



The conditions 



(J-v'h=b 2 

(f-vw=(br 



(14) 



(15) 
(16) 



and rj<b, rj' <b f may be used to determine rj and y f . 
Now for k greater than some number n replace 
j k by y\' for k even, by rj for k odd, and replace 
j k f by rj for k even, by -q' for k odd giving 



w k > 



(bb f ) 2 



and 



(J-v) 2 (j-vV 



W k '~pW k 2 



j- 2 W k ^2 = ^W k . 2 



(17) 



(18) 



for k sufficiently large. Thus the solution from (12) 
would be 



Wtt^iWk + W^ + lWn^ + Wn-i + Wn + W^T^ (19) 
k=0 1— P 

because 



*=i i- 



(20) 



An accurate solution for any value of bb'/f 2 can be 
obtained from eq (19) by an appropriate choice of n. 
It was found in the computations for the examples to 
follow that for bb / /f 2 =}U, n=l gave an accurate 
result, and for the extreme case bb / /f 2 = 1 /i an accurate 
solution was obtained for n= 3. 

Again when bb' 'If 2 is sufficiently small rj and t\ are 
negligible so that eq (19) reduces to 



W~ {Wo+Wi + Wo+W/} 



-m 



1 



(21) 



where w , w u w$ , and w x ' are given by eq (1), (2), 
(8), and (9) after setting f x and// equal to zero. 
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4. More Than Two Cylinders 

When more than two cylinders are moving in an 
infinite liquid the method of solution is completely 
analogous to that given in section 2. 

If n cylinders are moving in a liquid choose any 
one R as a reference and introduce an appropriate 
coordinate system. The complex potential of C R in 
the absence of the other cylinders will be that of a 
doublet w R , . The image of w R , is taken in each of 
the remaining n — 1 cylinders in the same manner as 
before. Each of these n— 1 image doublets will 
have an image in C R and in each of the other cylin- 
ders as well. The process is continued in this 
manner with each doublet introducing n—\ image 
doublets at each step. The same process is applied 
starting with each of the n—1 remaining cylinders. 
Again the potential of the system will be the sum 
of the potentials of all the doublets. An application 
of this method will be given in seel ion (>. 

5. Cylinder Near a Plane Boundary 

If b l = b, U'=U, a=ir—a the plane equidistant 
between the cylinder axes is a streamline and the 
solution represents the (low induced by a cylinder 
moving in the presence of a plane boundary. Figure 
2 shows the streamlines for a cylinder moving toward 
a wall, a=0. 

The equations can be simplified when the cylinder 
moves parallel to the wall, a=ir/2, with its circum- 
ference touching the wall, b/f =1/2. In this case the 
complete solution is 



w=iUbJ2 



r_ 



i 



f±i(k+l)U i 

[b k+l 



k+2 I 
k+l) 



(22) 



An approximate solution is obtained by taking n 
terms of the series and in the remaining terms re- 
placing k/(k+l) and (k+2) /(k+l) by 1. Thus 



n 1 r 
jfc=o (fc+l) 2 « 2 



1 



k+2 



+ iUb 



lb k+l b k+l 
2_ (V_A 



&- 1 



6 



§(fc+l) 2 



(23) 



The st reamlines are shown in figure 3a. The compu- 
tations were made 4 with yi=3. The streamlines for 
figure 3b were obtained by superposing a velocity 
— I r on the cylinder and the liquid. 

The streamlines for liquid (lowing around a cylin- 
der parallel to a plant 4 wall are shown in figure 4. 
Equation (19) was used with n=l, ol=t/2, 6//= 1/4. 




\\\\ \\\\\\\\\ 



Figure 



U=l,b=l. 



Stream lines for cylinder moving toward a plane 
boundary. 




Figure 3a. Streamlines for cylinder moving parallel to a plane 
boundary when the cylinder loaches the boundary. 
17=1, 6=1. 




Figure 3b. Streamlines for current flowing with general 
velocity U past a fixed cylinder against a plane boundary. 

Obtained from figure 3a by impressing a velocity — U on the fluid and the 
cylinder. 17=1, 6=1. 
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Figure 4- Streamlines for current flowing with general velocity 
U past a fixed cylinder near a plane boundary. 

17=1,6=1. 



6. Cylinder Near Two Plane Boundaries 

If four cylinders of equal radii have symmetric 
positions and velocities with respect to two perpen- 
dicular planes as indicated in figure 5, then each 
plane is a streamline and the solution represents the 
flow induced by a cylinder moving in the presence 
of two perpendicular plane boundaries. Consider- 
ing the case /= 4a, g=4a, a=0, an approximate solu- 
tion is obtained by the method of section 4 stopping 
after the first set of images. The streamlines are 
shown in figure 6. 







Q). 




Figure 5. Notation diagram for motion of a cylinder near two 
perpendicular plane boundaries. 

The dotted lines indicate the streamlines that can be considered as solid bound- 
aries. 




Figure 6. Streamlines for a cylinder moving toward one of two 
perpendicular plane boundaries. 
17=1,6=1. 



7. Cylinders Moving as a Rigid Body 

If the velocity, Ue ia , is the same for each cylinder, 
they move as a rigid body. The previous solution 
for a cylinder moving parallel to a plane boundary 
is one example of this case. Another example is 
given in figure 7 where b f = b=f/2, a ' = a=0. The 
solution simplifies to 



ptf (k+l)K 



1 



1 



z k_ ' z_k+2 

b k+1 b k+ 



?i 



(24) 



The approximate solution obtained by the method 
of section 5 is 



w~Ub±^-^(- 
k^o(k+l) 2 J z 



1 



1 



k 
b k+1 



+Ub 



z k+2 
b k+1^ 

12 fcrf(fc+l 



2_ (t\ ^ (-1)* ^ 
J 12 feUfc+1) 2 I 



(25). 



The computations were made with n=3. 

Another example is given in figures 8a and 8b 
where 6' = 6=f/4, a' = a=0. The solution was ob- 
tained using eq (19) with n=l. 
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Figure 7. 
17=1,6=1. 



Streamlines for current flowing with general velocity 
U past fixed cylinders. 




Figure 8a. 
U= 1,6=1. 



Streamlines for two cylinders moving as a rigid 
body. 




Figure 8b. Streamlines for current flowing with '[general 

velocity U past fixed cylinders. 
7=1,6=1. 
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suggestions of G. H. Keulegan and G. B. Schubauer. 



Washington, March 11, 1958. 
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